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We introduce a "hyperbicycle" ansatz for quantum codes which gives the hypergraph-product 
(generalized toric) codes by TiUich and Zemor and generaUzed bicycle codes by MacKay et al. as 
limiting cases. The construction allows for both the lower and the upper bounds on the minimum 
distance; they scale as a square root of the block length. Many of thus defined codes have finite rate 
and a limited-weight stabilizer generators, an analog of classical low-density parity check (LDPC) 
codes. Compared to the hypergraph-product codes, hyperbicycle codes generally have wider range of 
parameters; in particular, they can have higher rate while preserving the (estimated) error threshold. 



I. INTRODUCTION 

Quantum computing can become a reality only when 
combined with quantum error correction as quantum in- 
formation is intrinsically fragile due to inevitable cou- 
pling to the environment [IHl]. Quantum error correct- 
ing codes (QECCs) [5]-[7] offer protection for quantum 
information; however, often at a high cost in the number 
of auxiliary qubits and with technologically difficult re- 
quirements [8HT7L The optimization then must play an 
important role since by simplifying the syndrome mea- 
surements one can lower the overhead and perform er- 
ror correction in parallel. Quantum codes with limited 
stabilizer generator weights, a quantum analog of classi- 
cal low-density parity check (LDPC) codes [18, 19 seem 
to offer a viable solution here as the simple structure 
of stabilizer generators will require fewer ancillas, fewer 
rounds of syndrome measurements and parallelism. Fur- 
thermore, by analogy with classical LDPC codes, there 
might exist efficient algorithms for encoding and decod- 
ing [19]. 

However, quantum LDPC codes come at a higher price 
compared to their classical analogs due to stringent re- 
quirements on code parameters originating from commu- 
tativity of stabilizer generators. In fact, there are no 
known families of asymptotically good quantum LDPC 
codes, or any bounds suggesting their existence. It then 
becomes an intriguing question of the best asymptotic 
properties achievable with quantum LDPC codes. In 
such a setting explicit code designs become important, 
in particular, for establishing the lower bounds on the 
parameters: the number of encoded qubits K and mini- 
mum distance D for a given block length N (which also 
defines the code rate K/N)^ e.g., given the upper limit 
on the weight of stabilizer generators. 

The best-known quantum LDPC codes (that are also 
local) are Kitaev's toric codes and related surface codes 
with the minimum distance scaling as \fN [l0l[T71|20l[2T]. 
Existence of single-qubit-encoding LDPC codes with the 
distance scaling as a/TV log N has been proved in Ref. 23l: 
Tillich and Zemor proposed a finite-rate generalization 
of toric codes [23 . The construction relates a quantum 
code to a direct product of hypergraphs corresponding 
to two classical binary codes. Generally, thus obtained 
quantum LDPC codes have finite rates and the distances 



that scale as a square root of the block length. 

In one of the first studies of quantum LDPC codes 
MacKay et al. [19 constructed so called bicycle codes. 
Numerically, these codes exhibit good decoding proper- 
ties; however, the minimum distance of such codes is 
unknown. The quantum hypergraph-product codes [2^, 
on the other hand, are an example of LDPC codes with 
known parameters; however, decoding such codes may be 
difficult. The existence of a finite noise threshold, with 
and without syndrome measurement errors, has been re- 
cently established for limited-stabilizer-weight quantum 
hypergraph-product codes [24 , as well as for any such 
LDPC code family with the distance scaling as the square 
root of block length. These results, however, might 
not apply to the constructions of quantum LDPC codes 
based on finite geometries [25", ^26^ and to the construc- 
tions based on Cayley Graphs [27 due to the unbounded 
weight of stabilizer generators. 

In this work, we construct a large family of codes that 
in the limiting cases reduce to (generalized) bicycle and 
hypergraph-product codes, hence we call them hyperbicy- 
cle. We show that hyperbicycle codes contain new code 
families with finite rates and the distances that scale 
as a square root of the block length. In addition, the 
hyperbicycle construction can improve the rate of the 
hypergraph-product codes while preserving the estimated 
error threshold. 



II. PRELIMINARIES 

In this section, we define classical and quantum er- 
ror correcting codes. We also review some of the known 
LDPC code constructions. 



A. Classical error correcting codes 

A classical q-dLiy block error-correcting code {n^K^d)q 
is a set of K length-n strings over an alphabet with q 
symbols. Different strings represent K distinct messages 
which can be transmitted. The (Hamming) distance be- 
tween two strings is the number of positions where they 
differ. Distance d of the code C is the minimum distance 
between any two different strings from C. 
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In the case of linear codes, the elements of the al- 
phabet must form a Galois field F^; all strings form n- 
dimensional vector space F^. A linear error-correcting 
code [n^k^d]q is a /c-dimensional subspace of F^. The 
distance of a linear code is just the minimum weight of 
a non-zero vector in the code, where weight wgt(c) of a 
vector c is the number of non-zero elements. A basis of 
the code is formed by the rows of its generator matrix G. 
All vectors that are orthogonal to the code form the cor- 
responding (n — A:)-dimensional dual code, its generator 
matrix is the parity-check matrix H of the original code. 

For a binary code C[n, /c, d], the field is just F2 = {0, 1}. 
For a quaternary code C, the field is F4 = {0, 
with 



1, 



1, and uj ■. 



(1) 



For non-binary codes, there is also a distinct class of 
additive classical codes, defined as subsets of F^ closed 
under addition (in the binary case these are just linear 
codes). 

A code C is cyclic if inclusion (cq, Ci, . . . , c^-i) G C 
implies that (c^-i, cq, ci, . . . , Cn-2) ^ C. Codes that are 
both linear and cyclic are particularly simple: by map- 
ping vectors to polynomials in the natural way, c 
c{x) = C0+C1X+. . .-\-Cn-iX^~^ , it is possible to show that 
any such code consists of polynomials which are multiples 
of a single generator polynomial g{x)^ which must divide 
x'^ — 1 (using the algebra corresponding to the field Fg). 
The quotient defines the check polynomial h{x)^ 



h{x)g{x) = x^ 



(2) 



which is the canonical generator polynomial of the 
dual code. The degree of the generator polynomial is 
degg{x) = n — k. The corresponding generator matrix 
G can be chosen as the first k rows of the circulant ma- 
trix Cn formed by subsequent shifts of the vector that 
corresponds to g{x)^ explicitly: 
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B. Quantum stabilizer codes 



Binary quantum error correcting codes (QECCs) are 
defined on the complex Hilbert space ^ where is 
the complex Hilbert space of a single qubit a |0) + /3 |1) 
with G C and + = 1. Any operator acting 
on such an n-qubit state can be represented as a com- 
bination of Pauli operators which form the Pauli group 
of size 2^^+^ with the phase multiplier i^\ 



\I,X,Y,Z} 



m 



0,...,3, 



(4) 



where X, F, and Z are the usual Pauli matrices and 
/ is the identity matrix. It is customary to map the 
Pauli operators, up to a phase, to two binary strings, 
v,u G {0,1}^^ [28], 



U 



X-Z- ^(v,u). 



(5) 



^2 " '^n 



where = Xl^X^^ • • • and = Z^^^ J 
A product of two quantum operators corresponds to a 
sum (mod 2) of the corresponding pairs (v^,u^). 

An [[n,/c,(i]] stabilizer code Q is a 2^-dimensional 
subspace of the Hilbert space l-L®^ stabilized by an 
Abelian stabilizer group generated by the commut- 
ing Pauli operators (generators) Gi, . . . , Gn-ki i-e., ^ = 
(Gi, . . . , Gn-k) and -1 ^ y Explicitly, 



Q = m:S\^) 



(6) 



Each generator G ^ is mapped according to Eq. ([5| 
in order to obtain the binary generator matrix H = 
{Ax\Az) in which each row corresponds to a generator, 
with rows of Ax formed by v and rows of Az formed by 
u vectors. For generality, we assume that the matrix H 
may also contain unimportant linearly dependent rows 
which are added after the mapping has been done. The 
commutativity of stabilizer generators corresponds to the 
following condition on the binary matrices Ax and Az'^ 



AxA^^AzA^x =0 (mod 2). 



(7) 



A more narrow set of Calderbank-Shor-Steane (CSS) 
codes [30 contains codes whose stabilizer generators can 
be chosen to contain products of only Pauli X or Pauli Z 
operators. For these codes the stabilizer generator matrix 
can be chosen in the form: 



H 



Gx 






Gz 



(8) 



where the commutativity condition simplifies to 
GxGl = 0. 

The dimension of a quantum code is 



k = n — ranki^; 
for a CSS code this simplifies to 

k = n — rankGx — rankG^- 



(9) 



(10) 



The distance d of a quantum stabilizer code is given by 
the minimum weight of an operator U which commutes 
with all operators from the stabilizer but is not a part 
of the stabilizer, U y. In terms of the binary vector 
pairs (a, b), this is equivalent to a minimum weight of 
the bitwise OR(a, b) of all pairs satisfying the symplectic 
orthogonality condition. 



Axh^AzaL = 0, 



(11) 



which are not linear combinations of the rows of H. A 
code of distance d can detect any error of weight up to 
d — 1, and correct up to [d/2\. 
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In an equivalent representation, one can map any Pauli 
operator U in Eq. (|5|, to a quaternary vector over F4, e = 
u+cjv. A product of two quantum operators corresponds 
to a sum (mod 2) of the corresponding vectors. Two Pauli 
operators commute if and only if the trace inner product 
ei * 62 = ei • 62 + ei • 62 of the corresponding vectors is 
zero (which is equivalent to the symplectic orthogonality 
condition), where e = u + cJv. With this map, generators 
of a stabilizer group are mapped to rows of a generator 
G of an additive code over F4, with the condition that 
the trace inner product of any two rows vanishes [28 [see 
Eq. The vectors generated by rows of G correspond 
to stabilizer generators which act trivially on the code; 
these vectors form the degeneracy group and are omitted 
from the distance calculation. For CSS codes in Eq. (|8| 
the generator matrix is a direct sum G = Gx ® ujGz- In 
the following, we will use both, quaternary and binary, 
representations. 

A classical LDPC code is a code with a sparse parity- 
check matrix H. For a regular (j, /) quantum LDPC code, 
every column and every row of H have weights j and 
/ respectively, while for a (j, /)-limited quantum LDPC 
code these weights are limited from above by j and 

In the case of quantum LDPC codes, these properties 
apply to the generator matrix G whose rows correspond 
to the generators of the stabilizer group. 



C. Bicycle codes 

In one of the first studies of quantum LDPC codes, 
MacKay et. al. proposed a CSS code construction [19j 
which can be written in a block form as: 



G 



X 



G, 



{A A' 



(12) 



where A is a binary circulant matrix. Bicycle codes 
are obtained after some of the rows in Gx or Gz 
are deleted. Numerically, such codes show good error- 
correction capabilities [19| ,31^; however, the distance of 
such codes is unknown. 




FIG. 1: (Color online) Left: Two stabilizer generators 
(marked by arrows) and two pairs of anticommuting logical 
operators (marked by lines) of a [[450,98,5]] code in Eq. (13) 
formed by circulant matrices Hi = 7^2 corresponding to co- 
efficients of the polynomial h{x) = 1 + x + + (red - 
X operators, blue - Z operators, green - overlap of Z and 
X operators, dark and light gray - dual sublattices of phys- 
ical qubits). All other stabilizer generators are obtained by 
shifts over the same sublattice with periodic boundaries. In 
the shaded region, each gray square uniquely corresponds to 
a different logical operator, thus 98 encoded logical qubits. 
Right: same for the toric code [[450, 2, 15]]. 



(dimensions ri x ni) and I-L2 (dimensions r2 x 77,2), and Ei 
and Ei^ i = 1,2, are unit matrices of dimensions given 
by ri and n^, respectively. The matrices Gx and Gz^ 
respectively, have rir2 and nin2 rows (not all of the 
rows are necessarily linearly independent), and they both 
have N = r2ni + rin2 columns, which gives the block 
length of the quantum code. The commutativity condi- 
tion GxG^ = is obviously satisfied by Eq. ( |13| since 
the Kronecker product obeys {A^B){C^D) = AC^BD. 

The parameters [[A^, i^, D]] of thus constructed quan- 
tum code are determined by those of the four classical 
codes which use the matrices Hi, H2, , and I-L2 
the parity-check matrices. The corresponding parame- 
ters are introduced as 

Ch^ = [rii.ki.di], C^T = [ni,ki,di], i = 1,2, (14) 



D. Hypergraph-product codes 

Tillich and Zemor proposed a CSS construction which 
can be interpreted as a finite-rate generalization of toric 
codes [23 and allows for LDPC constructions. For such 
codes, the generator matrix is constructed from a product 
of two hypergraphs, each corresponding to a parity check 
matrix of a classical binary code. The associated CSS 
code can be recast in a matrix form with the generators 
given by [32 



Gx 
Gz-- 



= (^2' 

in! 



(13) 



Here each sublattice block is constructed as a Kronecker 
product (denoted with "(g)") of two binary matrices Hi 



where we use the convention [23 that the distance 
di{di) = CO if ki{ki) = 0. The matrices Hi are arbi- 
trary, and are allowed to have linearly-dependent rows 
and/or columns. As a result, both ki = rii — rank Hi and 
ki = Hi — rank Hi may be non-zero at the same time as 
the block length of the "transposed" code C^t is given 
by the number of rows of Hi, rii = ri. 

Specifically, for the hypergraph-product code ([l3|, we 
have N = r2ni + rin2, K = 2kik2 — kiS2 — k2Si with 

= ^i—^i-, ^ = 1,2 (Theorem 7 from Ref. [23 ), while the 
distance D satisfies the conditions D > min((ii, 6/2, <^i, <^2) 
(Theorem 9 from Ref. ^), and two upper bounds 
(Lemma 10 from Ref. [23]): if /ci > and /c2 > 0, then 
D < di; if A:2 > and ki > 0, then D < ^2- 

A full-rank parity check matrix Hi of a binary code 
with parameters C^i = [^1, ^1, di] (ri = ni — /ci, ki =0) 
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and I-L2 = define a quantum code with parame- 
ters [[(ni — ki)'^ + n^, /c^, (ii]][23 . Furthermore, a fam- 
ily of finite-rate (/i, 'L')-limited classical LDPC codes with 
asymptotically finite relative distance will correspond to 
a family of finite-rate {v^h -\- v)-limited quantum LDPC 
codes with the distance scaling as I) ex y/N. 



III. TWO-SUBLATTICE CODES 



specific for such codes [30^. Also, an original sparse code 
is mapped to a sparse code, with the same limit on the 
column weight, and row weight at most doubled. So, if 
one has a non-CSS code and wants to use one of the mea- 
surement techniques designed for such codes, this can be 
done by first constructing the corresponding CSS code. 

We use the reverse version of this mapping in Sec. |IVHl 
to construct the non-CSS versions of hyperbicycle codes. 



The commutativity condition for QECCs in Eq. ^ 
puts a strong limitation on suitable parity check matri- 
ces. The problem becomes even more difficult when the 
additional requirement of LDPC structure is imposed. 
In particular, this strongly limits possible counting ar- 
guments for establishing bounds on code parameters. In 
such a setting, constructions based on some ansatz be- 
come very useful. In the following, we study several CSS 
constructions based on two sublattices corresponding to 
the columns of the binary matrices 741(^2^) and Bi{A2): 



Gx = {Ai,B,), 



Gz = {Bj,Al) 



(15) 



where the matrices A^, 5^, i = 1,2, satisfy the condi- 
tion A1B2 -\- B1A2 = (we assume binary linear algebra 
throughout this paper). 



Two-sublattice CSS code from a generic 
stabilizer code 



A large number of two-sublattice CSS codes (15) can 



be obtained from regular stabilizer codes by the following 

Theorem 1. For any quantum stabilizer code [[N, D]] 
with the generator matrix 



H={A\B), 



(16) 



there is a reversible mapping to a two-sublattice quantum 
CSS code (1^ with Ai = = A, Bi = = B and 



the parameters [[2N,2K, D']], where D < D' < 2D. 

Proof. Explicitly, the generator matrices are 

Gx = {A,B), Gz = {B,A). (17) 

The dimension of the code simply follows from Eqs. ([9|, 
(10), given that rankGx = rankG^ = ranki^. Any 
binary vector e = (a|b) such that Ah + Ba. = maps 
to a pair of double-size vectors e^ = (b,a), = (a, b) 
which satisfy Gx^z = 0, Gz^x = 0; the corresponding 
weights obey the inequality wgt OR(a, b) < wgt(a, b) < 
2wgtOR(a, b), which ensures the conditions on the dis- 
tance. It is easy to check that the reverse mapping also 
works. □ 

Note that an original code that exceeds the generic 
quantum Gilbert- Varshamov (GV) bound [33 is mapped 
to a CSS code that exceeds the version of the GV bound 



B. Generalized bicycle codes 

Let us now start with two commuting square n by n 
binary matrices, AB -\- BA = 0. Then, we can satisfy the 
general two-sublattice ansatz ( 15 ) by taking Ai = A2 = 
A, Bi = B2 = B^ which gives 



Gx = {A,B), Gz={B^,A^ 



(18) 



In particular, the commutativity is guaranteed for cir- 
culant matrices A and 5, which corresponds to a gener- 
alization of the bicycle codes[19 , see Eq. (12). In this 
case, we map the linear combinations of rows in Gx to 
a classical length-n additive cyclic code over F4, where 
elements of the code are constructed from the gener- 
ator matrix G = ujA -\- B. If the circulant matrices 
A and B are generated by the polynomials fi{x) and 
/2(x), respectively, the space of the additive code corre- 
sponding to the stabilizer is generated by the polynomial 
g{x) = fi{x)uj + f2{x) modulo x^ — 1. That is, elements 
of the stabilizer are given by the coefficients of the poly- 
nomials h{x)g{x) mod(x"^ — 1), with arbitrary binary h{x). 

A canonical form of cyclic additive codes over F4 has 
been introduced in Ref.^3 where Theorem 14 states that 
any cyclic additive code can be represented via two gen- 
erators as {ujp{x) + q{x)^r{x)) mod{x'^ — 1) with p{x) = 
gcd[/i(x), — 1], r{x) = gcd[{x^ — l)f2{x)/p{x)^x'^ — 1] 
and degq{x) < degr(x) (gcd stands for the greatest com- 
mon divisor). In this case the code dimensionality is 
k = 2n — degp{x) — degr{x). The special case of cyclic 
additive codes with a single generator has been analyzed 
in Ref. [34] in which case the code dimensionality simpli- 
fies to k = n — degp(x), which formally corresponds to 
r{x) = x^ — 1. Note that, unlike the case of the usual 
quantum additive cyclic codes[28j|34], the mapping from 
Eq. ( p!8| ) works for any circulant matrices A and B; no 
additional commutativity condition is needed for the gen- 
erator polynomials p{x)^ q{x)^ and r{x). The parameters 
of thus obtained quantum codes are given by 



Theorem 2. The generalized bicycle codes in Eq. ([^ 
have the block length N = 2n, the number of encoded 
qubits K = 2degp{x) + 2degr(x) — 2n [K = 2degp{x) 
in the single generator case] and the distance exceeding 
or equal that of the classical code over F4 formed by code- 
words orthogonal to G with respect to the trace inner 
product. 

Proof. The code dimensionality immediately follows from 
the parameters of the canonical form of the code gener- 
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ated by g{x). The orthogonal code contains the quantum 
code, hence the distance estimate. □ 

Note that the distance estimate in Theorem [2] is tight 
only for pure codes since a possibility for degeneracy is 
not taken into consideration. 

Example 1. Suppose a cyclic linear code [n, /c, d] overW^ 
with a generator polynomial g{x) that divides — 1 gen- 
erates a code space . Then the parameters of the quan- 
tum CSS code in Eq. (T^ l are [[2n, 2n — 4deg g{x), > d]]. 
This construction is similar to non-CSS code construc- 
tion from linear cyclic codes in Ref 128], except that here 
the dual code does not have to he self- orthogonal. For a 
cyclic [30,25,4] code with g{x) = (1 + + ujx){l + 

X -\- ujx'^) we obtain a quantum code [[60,40,4]]. 



Several examples of such codes are given by the Haah's 
codes [35^, which give examples of local codes in 3D with- 
out string logical operators, and may lead to realizations 
of self-correcting quantum memories. Such codes are 
defined on two sublattices and have exactly the tensor- 
product structure discussed here. In Table [l| we list four 
codes presented in Ref. [35 . These codes are essentially 
constructed from a repetition code and it is straightfor- 
ward to generalize this construction to arbitrary cyclic 
binary codes by using the corresponding binary circulant 
matrix Hi. The commutativity of matrices A and B in 



Eq. ( 15 ) immediately follows. 



A non-CSS generalization of construction in Table [T| 
can be achieved by using symmetric circulant matrices 
[see non-CSS construction in Eq. (16)]. 



Example 2. A CSS family of odd-distance rotated toric 
codes fSEI is obtained for fi{x) = (1+x^^ +"^) and f2{x)_ 



18). 



x{l + x^^^ t = 1,2,... by construction in Eq. 
These codes have the parameters [[2t^ + 2(t + l)^,2,2t + 
1]]. Explicitly, [[10,2,3]], [[26,2,5]], [[50,2,7]], [[82,2,9]], 



The constructions in Theorems [T] and [5] [Eqs. ^17\ 
and (18) respectively] coincide for symmetric matrices, 
A = A^^, B = B^ . Then, starting with two palindromic 
polynomials fi{x) = x^^^^^^^"^ fi{l/x)^ i = 1,2, we can 



obtain non-CSS halved hyperbicycle codes in Eq. (16) by 
applying the reverse of Theorem [l] to the matrices A and 

B lag. 

Example 3. A non-CSS family of smallest odd-distance 
rotated toric codes 1341 is obtained for palindromic 
fi{x) = x\l + x^^'+i) and f2{x) = x'^\l + x^'"'^), 
t = 1,2,... by construction in Eq. (T^ . These codes 
have the parameters [[t^ + (t + 1)^, 1, 2t + 1]]. Explicitly, 
[[5,1,3]], [[13,1,5]], [[25,1,7]], [[41,1,9]], .... 

The codes in the last two examples exceed the lower 
bound in Theorem [2] due to degeneracy. 



IV. CSS AND NON-CSS HYPERBICYCLE 
CODES 

This section contains our most important results. We 
show that the families of hypergraph-product and gener- 
alized bicycle codes can be obtained as limiting cases of a 
larger family of hyperbicycle codes. The main advantage 
of this construction is that it gives a number of previ- 
ously unreported families of quantum codes with tight 
bounds on, or even explicitly known distance. This in- 
cludes many strongly degenerate LDPC codes with the 
distance much greater than the maximum weight of a sta- 
bilizer generator. In this section we discuss the construc- 
tion of such codes, their parameters, and give examples. 



A. CSS hyperbicycle codes: construction 



We define the hyperbicycle CSS codes as follows: 



n^- fx-LT^ 7(x) ^ ^ ^JM ^ 



(19) 



C. Tensor-product constructions and Haah's codes 

Further generalization of the bicycle-like construction 
in Eq. ([l8| can be achieved by combining tensor prod- 
ucts with commuting (e.g., circulant) matrices. The most 
general form of two-sublattice tensor-product codes has 
the form: 



...iu ^ii 



where Hf^ and Hf^ are matching, pairwise-commuting 
binary square matrices (i.e., Hf^H^^ + Hf^Hf^ = 0). For 
circulant matrices Hf^ and the commutativity is au- 
tomatically satisfied. 



Here we introduce two sets of binary matrices (dimen- 
sions ri X ni, z = 0, . . . , c — 1) and hi (dimensions r2 x n2, 
i = 0, ... ,c — 1); Ea^ El), Ea and ^5 are unit matrices 
of dimensions given by ri, r2, ni and 712, respectively. 
Matrices i^i^^) permutation matrices given by 



(x) 



a product of two permutation matrices, i.e. 
(l^^ = S^lJ) where {Ii)kj = ^j-k.imodc is a circulant 
permutation matrix, {S^)kj = (^j-/c,(/c-i)(x-i) mode and 
the positive integers c and x are coprime . A version of 
this construction for c = 2 and x = 1 has been previously 
reported by us in Ref. |^ 

The matrices Gx and G^, respectively, have crir2 and 
cnin2 rows (not all of the rows are linearly independent), 
and they both have 



N = c{rin2 + r2ni) 



(20) 
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Code 1. 


B^ni^ni^E + E^ni^Ui + ni^E^Hi 


Code 2. 


A^Hi^E^E + E^ni^E + E^Ui^ni + ni^E^ni+ni^ni^ni 
B^E^E^ni+ni^Hi^E + E^ni^ni + ni^E^ni 


Code 3. 


A ^ Ui ® E ® E + E ^Ui ® E + E ^Ui ^Ui +ni ® E ^Ui +ni ^Ui ^Ui 

B^Hi^E^E + E^ni^E + E^E^ni+ni^ni^E + E^ni^ni + ni^ni^ni 


Code 4. 


A^E^ni^Hi + ni^E^ni + E^ni^E 

B ^ Hi ® E ® E + E ® E ^Ui + E ^Ui ^Hi + Hi ^Ui ^Ui 



TABLE I: Tensor-product two-sublattice representation of Haah's codes corresponding to Eq. (15) where 1-Li is a circulant 
matrix corresponding to parity check polynomial p{x) = 1 + x of a repetition code, £^ is a unit matrix of the same dimensions 
with and the summation is mod 2. 



columns, which gives the block length of the quantum 
code. The commutativity condition GxG^ = is obvi- 
ously satisfied by Eq. (19) since the permutation matrices 
commute with each other. Note that for c = 1 and X = 1 



we recover the hypergraph-product codes in Eq. (13) and 
for Ti = rii = 1, i = 1,2, (i.e., and bi given by binary 
numbers) we recover the generalized bicycle code con- 
struction in Eq. (18). 



In order to characterize codes in Eq. (19) it is conve- 



nient to introduce the "tiled" binary matrices: 



-Hi 
Hi -- 



r(x) 
Ttx) 



^2 
^2 



For example, for c = 5 and x = 2 we have 



7(xj 



(21) 



( ai a2 as \ 

CL/i CL^ Cii CL2 CLs 

a2 as a4 a^ ai 

as ai a2 as a^ 

\ as a^ a^ ai a2 ] 



(22) 



note that the subsequent block rows are shifted by x = 2 
positions. 

For the following discussion it is also useful to define 
auxiliary binary matrices: 

^1 = Y.ih ® au nl = Eibi ® li, (23) 
where in terms of these matrices we can write: 



(24) 



With this notation, it is clear that the generator matri- 
ces (19) correspond to the hypergraph generators (13), 



except that the identity matrices Ei^ Ei are now reduced 
in size by a factor of 1/c. The block length of the original 
hypergraph code with the present binary matrices (|2T|) is 



orig 



\r1n2 +r2ni). 



(25) 



In this sense, one can view the codes defined by Eq. (19) 
with c > 1 as reduced hypergraph codes. 



B. CSS hyperbicycle codes: dimension 



Just as for the hypergraph codes, the parameters of 
classical codes C-^. and with parity check matrices 

1-Li and Hi, respectively, contain information about the 
parameters of the quantum code in Eq. (19). We denote 



the distances of these binary codes as di and d^, and their 
dimensions as ki and /c^, i = 1,2. 

Regardless of the choice of the matrices a^, 6^, the codes 
C^i and Cy^ are quasicyclic, with the cycle of length 
equal the dimension c of the cyclic permutation matri- 
ces I^i {Ixi)i ^ = l?^. Indeed, the corresponding block 
shifts merely lead to permutations of rows of the check 
matrices Hi^ Hi [see Eq. (22)]. In order to define the 
dimension of the corresponding hyperbicycle codes with 
generators (19), we first classify the vectors in C^i and 
Cy^ with respect to this circulant symmetry. 

We start with the case of a binary cyclic code with 
block length c, with the generator polynomial g{x)^ which 
divides — 1. [Polynomial algebra in this section is done 
modulo 2.] Any codeword corresponds to a polynomial 
w{x) which contains g{x) as a factor, and, therefore, ev- 
ery cyclotomic root of g{x) is also a root of w{x). How- 
ever, the particular polynomial w{x) = g{x)f{x) may 
also contain other factors of — 1 and thus have sym- 
metry different from that of g{x). We can define a linear 
space C^^^ of length-c vectors corresponding to w{x) with 
the exact symmetry of g{x) where p(x) = {x^ — l)/g{x) by 
defining the equivalence wi{x) = W2{x) for a given g{x) 
as /i = /2 modp^{x) for all p'{x) such that p'{x) ^ p{x) is 
a factor ofp{x). The same equivalence can be also defined 
modulo greatest common divisor (gcd) of all such poly- 
nomials p'{x). In terms of the corresponding check poly- 
nomial p(x), the dimension k^^^ of thus defined space C^^^ 
is zero unless p{x) is a non-zero power of an irreducible 

polynomial Pa^^)-, in which case k^^^ = degpa{x). 

For the quasicyclic code with the first check ma- 
trix in Eq. (21), the vector w is in the symmetry class of 
p(x), where p(x) divides — 1, if w satisfies the condi- 
tion [p{Ii) (g) 

^a] w — and is not a member of such a 
symmetry class of any factor of p{x). For each polyno- 
mial p{x)^ the l.h.s. in these equations is a sum of cyclic 
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shifts of the vector w corresponding to each non-zero co- 
efficient ofp{x). We denote the dimension of the subcode 
of with ah vectors in the symmetry class of p{x) as 



The symmetry imphes that k\^^ must contain the 
dimension k^^ introduced in the previous paragraph as 
a factor, and, in particular, kl^ must be zero whenever 



b) 



is zero. A convenient basis of C^^ can be constructed 



using the following 



Hi 



Lemma 1. Any vector of the subcode can he chosen 
in the form 



w 



i=0 



(26) 



where the vector g corresponds to the generating poly- 
nomial g{x) = {x^ — l)/p{x); the vectors {Is (8) Ea)^^, 
< 5 < /cq^\ are linearly independent. 

Proof. Any vector of the subcode can be expanded 
in the form uj\^^ (8)6^, where are all distinct weight-one 
vectors, and uj\^^ are vectors from C^^\ Generally, any 
vector u; G C^^^ can be written as a sum 



of shifts of the 

vector g, Xlslo ^ ^^§5 obtain Eq. (26) by rearranging 



the summations. Linear independence follows from the 
symmetry of the vectors cj-^^ ^ C^^\ □ 

Note that, in addition to the symmetric vectors, the 
code Cq-i-^ may contain vectors with no special symmetry 
with respect to the discussed block shifts. We will for- 
mally assign these to the check polynomial p{x) = — 1, 

and define k^^ = 1. 

For the vectors of the code Cq-i^ with the second check 
matrix in Eq. ( [21] ), the condition to be in the sym- 
metry class of p(x) reads [£^5 (g) p(/i)]w = while 
[E}) (g)p'(/i)]w 7^ for all p'{x) 7^ p{x) that divide p{x). 

We denote the dimension of the corresponding subcode 

(v) 

as k2 ' The same classification can be done for codes 
with the transposed check matrices; the correspond- 
ing dimensions are k^j^\ i = 1, 2. 

The introduced symmetry classification is in the heart 
of the following 

Lemma 2. A vector v that belongs to both Ce^^Ui 
and 0^2^ Ea 'rnust he in the same symmetry class p{x) 
with respect to both codes Hi and I-L2, including the no- 
symmetry case p{x) = x^ — 1. Any such vector can he 
generally expanded in terms of 



Proof. The parameter x does not enter this discussion 
since it corresponds to permutations of rows in matri- 
ces and That the symmetry must be the same 
becomes evident if we write the most general expansion 



(28) 



where are length-c vectors and e-^ ^ are distinct 
weight- 1 vectors. Indeed, the condition to be in the 
symmetry class of p{x) is the same for both codes: 
[El) (g) p{Ii) (8) Ea]v = for p{x) itself but not for an y of 
its factors; thus 7^^ must be in C^^^. The expansion (^27 
follows from Lemma [T] □ 

We can now count linearly-independent rows in the 
generator matrices: 

Lemma 3. The numbers of linearly independent rows in 
matrices (l^ l are 

rankGx = rir2C -^k['''%'''^ /k^J''\ 
I 

rankGz = nmsc - ^ i^^^^^i^^^^V4'''\ (29) 



where pi{x) are all binary factors of x^ — 1 such that 
k^Q^^ 7^ 0; including x^ — 1 itself. 



Proof. We first count linearly-dependent rows in Gz- 
Notice that the equations • (^5 Hi) = and 
^ ■ {H2 ^ Ea) = are both satisfied for v in Eq. (27) 



V 



[Lemma [2]. Each pair (a,/3) generates k^^^ linearly- 
independent vectors, same as each of them generates for 
the corresponding subcodes i = 1,2, respectively. 

Thus there are exactly k[^^ k2^^ / k^^^ linearly-independent 
vectors corresponding to every p{x) with non-empty C^p\ 
Such vectors have to be complemented with the pairs of 
vectors of no symmetry (if any) which formally corre- 
spond top{x) = x^ — 1 and /cq^^ = 1. According to Lemma 
|2] these are all possible solutions, which gives rankG^ in 
Eq. (29). We obtain rankGx by substituting the pa- 



rameters of the codes with the parity check matrices Hi, 
H2- □ 

We finally obtain 

Theorem 3. A quantum CSS code with generators 
encodes 



(27) 



where ^ • 13% ^ lig ^ C^l and ^ • /^g G cj^^^ and g 
corresponds to the polynomial g{x) = {x^ — l)/p{x). 



K = 2 ^ i^f '^4^^^ - kiS2 - k2Si (30) 



qubits, where pi{x) are all binary factors of x^ — 1 such 

that /cq^'^ 7^ 0; including x^ — 1 itself and Si = rii — ri, 
i = 1,2. 
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Proof. The number of encoded qubits K can be deduced 
from Lemma [3] using the relation 



k 



ip) 



Sjk, 



iP) 



1,2. 



(31) 



The latter follows from the fact that the rank of a ma- 
trix does not change under transposition (and also un- 
der permutations of rows and columns, e.g., as needed to 
transform 1-Li into 1-Li). Specifically, restricting the action 
of matrices l-Li and l-Li to subspace C^'p\ we obtain re- 
duced mutually transposed matrices of dimensions given 



by Vik^Q^ X n^/cQ^^ and nik^^^ x r^A:, 



■Sp) 



{p) 



gives Eq. (31 ). 



which immediately 

□ 



By construction, any k\ may only be non-zero if the 



corresponding /c^ > 0, i = 1,2. Then, Eq. (30) gives 



Consequence 4. A quantum CSS code with genera- 
tors (1^ can only have K > if at least one of the binary 
codes with the parity check matrices (21) is non-empty. 



The upper bound on the distance of the code with 



generators (19) is formulated in terms of thus introduced 



subset-distances d) 



ip) 



d^'\i = 1,2: 



Theorem 6. For every p{x), a binary factor of — 1 



such that k^^ > and k)^^ > 0, the minimum dis 



Zip) 



tance D of the code with generators /[ISy satisfies D < 
min(4''\5^/^) 



Similarly, when k2^^ > and k^^^ > 0, 



we have D < mm{d2^\ d^^"^ ) . 

Proof. Given k^^ > 0, consider vector u = (e (g) c,0). 



where c G cj^^^ and wgt(e) = 1. As long as k2^^ > 0, we 
can always select such e that u is not a linear combination 
or rows of Gz, which would indicate that D < wgt(c). 

Indeed, by construction, vector c can be written in the 
form ( 26 ) ; let us pick a bit s which is not identically zero 



in all ai and construct a vector [cf. Eq. (p^ 



(35) 



C. CSS hyperbicycle codes: general distance 
bounds 



Theorem 5. The minimum distance of the code with 
generators (19) satisfies the lower bound 



Taking all r2 different vectors e and all k^^ linearly- 
independent translations [Lemma l], we obtain the vector 
space [as indicated in Eq. (|35| with a brace] isomorphic to 

that on which the subcode C^~} operates. On the other 

rt2 



D > [d/c\, d = mm{di,d2,di,d2). (32) hand, there are only r2k^^^ - k^^^ linearly-independent 



:ip) 



Proof. Consider a vector u such that Gx • u = 0. We 



combinations of rows of the matrix I-L2 restricted to the 

Zip) 



construct a reduced quantum code in the form (IT9|), with subspace C^p\ Since A:/ > 0, at least one of vectors 



the same c, by keeping only those columns of the matrices 
a^, hi that are involved in the product Gx • u. This way, 
for every non-zero bit of u, one of the reduced matrices 
l-L'2 [see Eq. ( [21] )] may get c columns, so that these 
matrices have no more than cwgt(u) columns. If we take 
wgt(u) < \_d/c\ , according to Consequence [4j the reduced 
code encodes no qubits, thus the corresponding reduced 
u', -u' = 0, has to be a linear combination of the rows 
of G'^. The rows of G'^; are a subset of those of G^, with 
some all-zero columns removed; thus the full vector u is 
also a linear combination of the rows of Gz- Similarly, a 
vector V such that Gz • v = and wgt(v) < \d/c\, is a 
linear combination of rows of Gx- D 

Let us introduce the minimum distances d^^^ corre- 
sponding to the subset of the vectors of the code C-^. 
which contain one of the vectors with the exact symme- 
try of p{x), 



d^^ = min{wgt(a + b)|0 7^ a G C^^/, b G C^A C'^!}- 



ip) 



AP)^ 



Evidently, thus introduced distances satisfy 

^iP) rl. r^\r. rliP^) 



dY^ >di, mind]'''^ =di, z = l,2; 



(33) 



(34) 



the minimum is taken over all pi{x) as in Theorem Is] 
We will also introduce the distances d]^ corresponding 
to the matrices Hi, i = 1,2. 



uip),s ig 

linearly independent of the rows of the matrix 
I-L2 restricted to the subspace C^p\ 

Now, we can construct such a vector u for every c from 
the set in Eq. (34), which proves D < d^^\ The other 
bounds in the Theorem can be obtained from this one by 
considering isomorphic codes [e.g., interchanging I-L2 and 
Hi, and also Hi and H2]. □ 

The meaning of the condition on p{x) in Theorem [6] 
can be elucidated if we rewrite the number of encoded 



qubits (|30|) with the help of identity (31), 



K iLbOiLbO iT^iPi) I h^P'H^P'^ Ih^' 

j\ — 2_^i^i 1^2 / '^0 ^ 2 ^1 /^o 

I 



ipl) 



(36) 



Obviously, every term in Eq. ( 36 ) giving a non-zero con- 



tribution to K, also gives an upper bound on the mini- 
mum distance of the quantum code. 



D. Codes with finite rate and distance scaling as 
square root of block length 

Here we show that the family of hyperbicycle codes 
contain (v, h + 'u)-limited LDPC codes with the distance 
D (X >/]V that are distinct from the hypergraph product 
codes. Let us start with a random {h, v)-regular parity 
check matrix of a classical LDPC code, where h < v. 
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By removing linearly dependent rows, we can form full- 
rank (/i, 'u)-limited parity check matrix ai that, along 



with 6. 



(j 7^ 1 and x arbitrary, in case when 



j = I and X = 1 we recover the hypergraph-product 
codes), we use in Eq. (19) in order to construct the hy- 



perbicycle code where only one term in each summation 
in Eq. (19) is taken. The rate of the classical code defined 



by the parity check matrix ai is bounded from below, i.e. 

= kc/ric > l — h/v. With high probability at large nc, 
the classical code will also have the relative distance in 
excess of some finite number ^dSS]. If the classical LDPC 
code defined by ai has parameters then, ac- 

cording to Theorem [3) [6] and [5) the quantum code will 
have parameters [[c{ni — /ci)^ + cn^, c/c^, > di/c]]. It fol- 
lows that a finite rate (/i, i')-limited classical LDPC code 
(defined by the parity check matrix ai) with finite rela- 
tive distance (we expect the subset relative distance in 
( [34| ) to be finite as well) will correspond to a finite rate 
(vTh + 'L')-limited quantum LDPC code with the distance 
D (X V^. 



E. Codes with repeated codewords 

In some cases the distance of the hyperbicycle codes is 
larger than the lower bound in Theorem [5] In this section 
we consider the special case of square matrices a^, bi {vi = 
ni), with the additional restriction that the codes C-^^, 
are non-empty {ki = ki > 0) and contain only fully- 
symmetric vectors in the symmetry class of p{x) = l-\-x. 
The results we proved so far give the parameters of such 
codes summarized by (see also Theorem 3 in Ref. ^3^ ) 



Consequence 7. Suppose ai and bi in Eq. ll21y are such 
that 



ki > and Vi = n^. Then the CSS code 
with generators JT^ has the block length N = 2cnin2, 
encodes K = 2kik2 qubits, and has the minimum distance 
D limited by \d/c\ < D < d, d = m.m{di , (i2 , g^i , ^^2 ) • 

Proof. By assumption, all vectors in the codes C-^., i = 
1,2, are in the symmetry class of p{x) = 1 + x, which 



corresponds to k^ 
the form 



1 and block-symmetric vectors in 



wi = g (g) a, W2 = /3 ( 



(37) 



respectively, with g = (1, . . . , 1) [see Eq. (26)]. The num- 
ber of encoded qubits K immediately follows from The- 
orem [sj the block length N from Eq. ( pQ|) , and the lower 
bound on the distance from Theorem [5J Furthermore, 
with all vectors in the binary codes having the same sym- 
metry, the upper bound in Theorem [6] is just D < d. □ 

At this point we notice that the proof of the lower 
bound ld/c\ on the distance in Theorem [s] implies that 
there may be uncorrectable errors of the form "^gi/Ss ^ 
gs /35, a'^ gs cxs), where ah g^ have wgt(g5) = 1. 
On the other hand, if we were to consider only fully- 
symmetric vectors, with gg = g = (1,...,1), the factor 
of 1/c would be unnecessary. We formulate this result as 



Lemma 4. A symmetric vector u = (wi,W2); = 
/3* (g)g(g)a* with g = (1, . . . , 1); z = 1, 2^ that satisfies 
Gx^ = and is linearly independent from the rows of 
Gz, has sublattice weights wgt(wi) either zero or > d. 

Let us first consider the case c = 2 (then x must be 
equal to 1); we previously formulated the sufficient con- 
ditions to increased lower distance bound as Theorem 3 
in Ref. [32^ which was given without a proof. 



Theorem 8. Suppose c = 2, ai and bi in Eq. (21) are 



such that k^^^^^ = ki > 0, ri = rii and binary codes with 
generator matrices X] XI ctnd ^ bj have dis- 

tances at least 2. Then the CSS quantum code with gen- 
erators Eq. (19) has parameters [[4nin2, 2/ci/c2, where 
d = mm{di^d2,di^d2)- 

Proof. In addition to what is stated in Consequence [7[ 
we only need to prove that d is also the lower bound on 
the distance. To this end, notice that any vector u such 
that Gxu = can be decomposed as the sum of an "ac- 
tual" solution plus degeneracy, u^^+^^ + l^Gz^ where 
= (wi,W2) is a block-symmetric vector satisfy- 
ing the conditions of Lemma [4] and linearly-independent 
from the rows of Gz- This decomposition can be verified 
by comparing K with the number of linearly-independent 
solutions in the form (35), as well as those on the other 



sublattice. First, let us assume wgt(wi) > and there- 
fore wgt(wi) > d. We can rewrite the corresponding 
decomposition as wi = /3s (g) g (8) e J , where g = (1, 1), 
each e\ has length ni and wgt(ej) = 1, with the non- 
zero element in the position 5; there must be at least d/2 
non-zero vectors j3s. The full solution including the de- 
generacy can be formally written as w'^^ (g) ej, where 



<s = /3s ® g + 7X1, 0) + 7" <^ (1, 0), 



(38) 



where the sum of the last two vectors is a linear combina- 
tion of rows of I-L2 ' The key to the proof is the observation 
that 7^ + is a linear combination of rows of a^ + , 
and therefore is in the binary code generated by a^; 
by condition the corresponding weight is either zero or 
> 2. Without limiting generality we can drop the case 
7^ = 'jg which corresponds to a symmetric vector 
and can be included as a part of u'^^^^^. We are left with 
the trivial 7^ = 7^' = 0, in which case w'^^ = /3s (g) (1,1) 
remains unchanged; otherwise 7^ 7^ 7^^, in which case the 
weight of the modified w^^ can be lower bounded by that 
of the sum of the components corresponding to (1,0) and 
(0,1), 



wgt(w;j >wgt(7;+7;') >2; 



(39) 



with at least d/2 such terms the total weight is d or 
greater. The same arguments can be repeated in the case 
wgt(w2) 7^ 0, as well as for the space orthogonal to Gx- 
Overall, this proves the lower bound D > d; combined 
with the upper bound we get D = d. □ 
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Theorem 9. Suppose c is even, ai and hi in Eq. \21y 

are such that k^^^^^ = ki, ri = ni and binary codes with 
generator matrices ^ai, , and bj have dis- 
tance at least 2. Then the quantum code in Eq. jT^ has 
parameters [[2nin2C, 2/ci/c2, 1^]]; where {2/c)d < D < d 
and d = min((ii, (i2, di^ d2)- 

Proof. The proof is similar to the proof of Theorem [Sj 
except that now vectors w'^^ are defined by the analog 
of Eq. (38) which has g = (1, . . . , 1) with c components 

and more terms with 7s in the r.h.s., j = 1, . . . , c. We 
need to show that a non-zero w'^^ has wgt(w'^g) > 2, 
which ensures that the minimum distance of the code is 
at least 2d/c. 

With c > 2 and even, after the summation over all 
possible shifts of the vector w'^^ with respect to the block 
structure the symmetric term disappears, and we obtain 
the inequality cwgt(w'^g) > cwgt(7^ + 7^ + . . . + 7^). 
The sum in the r.h.s. is a linear combination of rows of 
^af; by assumption, it's weight is either > 2 or zero. 
The only non-trivial situation corresponds to the latter 
case with some 7^^ 7^ 0. For the sum to be zero, either 
there is an even number m of identical vectors 7^^ = 
= . . . = 7^^ ^ with m < c and all indices different 
[this situation results in wgt(w^g) > (c — m) > 2 since 
both m and c are even and l3s 7^ 0], or there are at least 
two pairs of unequal vectors 7^^ ^ 7^^ and 7^^ 7^ 7^^, 
with 7f2 ^ 7^4^ which also gives wgt(w'^g) > 2. □ 

In order to obtain codes with repeated structure (see 
Fig. |2|, one can start with two cyclic LDPC codes with 
block lengths n^, i = 1,2, and the check polynomi- 
als hi{x) that divide x^' — 1. The polynomials hi{x) 
will also divide x^'^' — 1, thus the corresponding circu- 
lant parity-check matrix Hi of dimensions cui x cui will 
lead to a code with repeated structure satisfying The- 
orem [9] since the corresponding generator polynomial is 

gi{x) = (x(^-l)^^ + ^(c-2)n, + . . . + 1) (^n, ^ l)/hi{x), 
i = 1,2. 

Example 4. Suppose we use the polynomial h{x) cor- 
responding to the shortened Reed-Muller cyclic code with 
parameters [2^ — l,m + 1,2^"-^ — 1] in order to con- 
struct circulant matrices Hi = I-L2 of dimensions 2(2^ — 
1) X 2(2^ — 1). According to Theorem a code in 
Eq. (1^ with = 2 and x = 1 will have parameters 
[[4(2^ - 1)2, 2(m + 1)2,2(2^-1 - 1)]]. This family leads 
to weight limited LDPC codes and up to m = 11 there 
is always a choice of polynomial h{x) of weight 4 which 
leads to quantum LDPC code with stabilizer generators 
of weight 8. 

Example 5. Given two ''small" cyclic codes [rii^ki^di] 
with check polynomials hi{x), i = 1^2, we can construct a 
0=1 hypergraph-product quantum code with the parame- 
ters [[2nin2, 2/ci/c2, d]], d = min((ii, d2), a repeated even-o 
code with the parameters [[2onin2^2kik2^D]], 2d < D < 
od, or a hypergraph-product code [[2o^nin2^ 2/ci/c2, do]] us- 
ing the ''large" cyclic codes with the same check polyno- 
mials and the block lengths oui. 



Note that in this Example the code rate goes down 
compared to the hypergraph-product code constructed 
from the "small" cyclic codes and goes up compared 
to the hypergraph-product code constructed from the 
"large" cyclic codes. 



11: 








FIG. 2: (Color online) Same as Fig. [T] for the [[294,18,4 < 
D < 12]] code in Eq. (19) formed by circulant matrices 
Til = l-i2 corresponding to coefficients of the polynomial 
h{x) = 1 + X + with c = 3 and x — 1- Two stabilizer 
generators are marked by red and blue arrows, respectively, 
and two anticommuting logical operators are marked by red 
and blue lines, respectively. All other stabilizer generators 
are obtained by shifts over the same sublattice with periodic- 
ity in the horizontal direction and shifted periodicity (shown 
by arrows) in the vertical direction. In the shaded region, 
each gray square uniquely corresponds to a different logical 
operator, thus 18 encoded logical qubits. One can observe 
the tripling of the logical operators, thus the overlap (green 
square) is also repeated three times. 



Planar qubit layout of hyperbicycle codes and 
encoding 



The stabilizer generators corresponding to Eq. ( 19 ) can 



be graphically represented on two rectangular regions 
corresponding to two sublattices. In case, when matri- 
ces Hi and I-L2 are square, the rectangular regions of 
sublattices have the same dimensions and can be drawn 
together with parameters and x corresponding to the 
number of square blocks and boundary shift, respectively, 
see, e.g.. Fig. |3] Furthermore, in some cases, we can 
represent logical operators by line-like operators with a 
possibility of using this layout for encoding. 

We start by considering the case o=l and x = 1 cor- 
responding to the hypergraph-product codes. The sta- 



bilizer generators for the quantum code in Eq. (13) can 



be graphically represented by two (dotted) lines living on 
different sublattices with the dots (red and blue squares 
in Fig. 1] marked by arrows) placed in the positions cor- 
responding to Is in rows of the binary matrices Hi, Hoq 
Hi = Hi and H2 = Hj - For cyclic codes, e.g., in Fig. Ill 
the relative position of dots stays the same and we can 
translate each stabilizer generator over the corresponding 
sublattice. In general, the form of stabilizer generators 
is position dependent and the peculiar two-line structure 
(see Fig.JT]) ensures commutativity. The logical opera- 
tors Xj, Zj^ j = 1, /c can be chosen among the rows of 



the matrices A'l 



{H^^^Ei.O), X2 = {0,E2^Hj^) 
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FIG. 3: (Color online) Upper plot: visualization of a [[90, 2, 9]] 
toric (hyperbicycle) code in Eq. ( 19 ) formed by circulant ma- 
trices 7^1 and 7^2 corresponding to coefficients of the polyno- 
mial h(x) — 1 + X, c — 5 and x = 3. The boundaries are 
periodic in the horizontal direction and shifted (as shown by 
arrows) in the horizontal direction by x = 3 blocks. Lower 
plot: general block construction leading to rotated periodic 
boundaries of hyperbicycle codes for c = 5 blocks and for the 
shift X = 3. This corresponds to t = 1 case of the infinite se- 
ries of block constructions with c = t^ + (t+l)^ and x = 2t+l, 
e.g., in case of toric code we get [[2n^c, 2, nx]] for any integer 
n > 1. 



bilizer generators is position dependent in case of non- 
cyclic codes. In general, codes with c > 1 have compli- 
cated structure of logical operators. Nevertheless, in a 
specific case of CSS codes when c is odd and k^^^^^ = ki 
(see Theorenij3|, we can recover the form of logical op- 
erators Xj, Zj, j = l,...,/c obtained in c = 1 case 
where the operators can be chosen among the rows of 
the matrices A^i = (H^ ^i,0), ^2 = (0,^2 ^ ^i) 
and Zi = {E2 (g) H^,0), Z2 = (0,H^ (g) ^1). The only 
difference is that the logical operators are now repeated 
c times which can lead to codes with increased distance 
(see Fig. [2|. 

Example 6. A CSS [[900,50,14]] hyperbicycle code is 
obtained with circulant Hi = Hi corresponding to the 



polynomial h{x) = (1 
X = 1; and hi = ai. 



-\- X -\- X -\- X ), Hi = 15; c = 2, 



Note that one-to-one correspondence between a set of 
physical qubits (shaded region in Figs, [l] and [2]) and 
logical qubits can be used for encoding. 



and Zi = {E2 (g) H]^,0), Z2 = (0,H^ g) £^1) where the 
index corresponds to the sublattice number on which the 
logical operator lives, ± stands for the orthogonal space 
mod 2 and matrices H^, and ^2^^ are in a 

row echelon form. By row and column permutations on 
matrices Hi, I-L2 it is convenient to reduce matrices Tii, 
^2": and 1-L2^ to the form with an identity matrix 

on the right. In such a case, the logical operators can be 
represented by vertical and horizontal (dotted) lines that 
have only one non-zero element in the region of the size 
ki X /c2 for the first sublattice and of the size ki x k2 for 
the second sublattice (shaded region in Fig. [T]) result- 
ing in k = kik2 -\- kik2 logical qubits. Thus, for such a 
representation, each physical qubit in the region of size 
kik2 + kik2 (shaded region in Fig. [l]) overlaps with only 
one logical qubit and can be used for encoding. Note that 
in general the two sublattices cannot be drawn together 
as they will have different dimensions for non-square ma- 
trices 1-Li and H2. In such a case, the sublattices can be 
represented by two different rectangular regions and the 
stabilizer generators have one line per sublattice. 



The hyperbicycle construction in Eq. (19) for arbi- 
trary c and X has a block structure of several rectan- 
gular regions stitched together with one of the periodic 
boundaries being shifted by x blocks (see Fig. [3|. The 
stabilizer generators can be graphically represented by 
two (dotted) lines with the dots (red and blue squares 
in Fig. [4| placed in the positions corresponding to Is 
in rows of the binary matrices Hi, H2, Hi and H2. 
For cyclic codes, e.g., in Fig. [4j the relative position 
of dots stays the same and we can translate the stabi- 
lizer generator with (shifted) periodic boundaries. Just 
like for the hypergraph product codes, the form of sta- 
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FIG. 4: (Color online) Same as Fig. [1] Left: X and Z sta- 
bilizer generators for the CSS hyperbicycle code [[900, 50, 14]] 
formed by circulant matrices corresponding to coefficients of 
a polynomial h{x) = 1 + x + + and c = 2, x = 1- 
Right: a single stabilizer generator of a [[289,81,5]] non- 
CSS hyperbicycle code in Eq. ( 16 ) formed by circulant ma- 



trices corresponding to coefficients of a polynomial h{x) — 
1 + X + + + + and c — 1. The division into two 
sublattices is impossible and all other stabilizer generators are 
obtained by shifts over the light and dark gray qubits with 
periodicity in the vertical direction and shifted periodicity 
(shown by arrows) in the horizontal direction. 
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G. Codes from two circulant matrices 



The hyperbicycle construction in Eq. ( 19 ) can employ 



the known families of cyclic codes when matrices and 
I-L2 (after additional permutations that change the or- 
der in the Kronecker product) in Eq. (23) are circulant. 



Note that any circulant matrix will have the block form 
of Eq. (23). As was mentioned in the previous section, 
for circulant matrices Hi and I-L2 the stabilizer genera- 
tors are translationally invariant with (shifted) periodic 
boundaries. 

The choice of x 7^ 1 can lead to codes with increased 
distance. This can be best seen on the example with the 
toric code (Fig. |3| where by rearranging the surface of 
the code we can bring it into a new layout with proper 
periodic but rotated boundaries. Then the Manhattan 
distance (defined on blocks, e.g., 3 in Fig.|3| between the 
boundaries will actually determine the distance of the 
code. The largest distance can be expected for squares 
thus by defining the boundary angle as {t^t -\- h) we ar- 
rive at codes with c = -h (t -h h)'^^ x = 2t + and 
the Manhattan distance equal to x- Compared to the 
general distance bound in Theorem [5] for toric codes we 
achieve the distance: D = x^/c. As the following exam- 
ples confirm, numerically we see that x > 1 can produce 
codes exceeding the distance bound in Theorem [5l often 
saturating the upper distance bound in Theoreni|6j 

Example 7. A CSS family of rotated toric codes is ob- 
tained when Til corresponds to the polynomial h{x) = 
(1 + x) (for Ul we use hi = ai), c = t^ ^ (t ^ 
X = 2t -\- 1, t = 1,2,.... By construction in Eq. (l^ 
we obtain codes with parameters [[2n^c, 2, nx]] . Explic- 
itly for n = 2 we obtain [[40,2,6]], [[104,2,10]]..., and 
forn = 3 [[90,2,9]], [[234,2,15]] .... 

Example 8. A [[90,8,8]] CSS hyperbicycle code is ob- 
tained when Hi corresponds to the classical cyclic code 
[15,4, 8] with the generator polynomial g{x) = (1 + + 
x^) (for I-L2 we use hi = ai), c = 5 and X = 3. 

Example 9. A [[90,10,7]] CSS hyperbicycle code is ob- 
tained when Hi corresponds to the classical cyclic code 
[15,5,7] with the check polynomial h{x) = {l-\-x-\-x^-\-x^) 
(for I-L2 we use hi = ai), c = 5 and X = 3. 

Example 10. A [[126,8,10]] CSS hyperbicycle code is 
obtained when Hi corresponds to the classical cyclic code 
[21, 5, 10] with the check polynomial h{x) = (1 + x + x^) 
(for Hi we use hi = ai), c = 7 and x = 3. Same con- 
struction with X = 1 results in the code [[126, 14,6]]. 

Example 11. Same construction starting with the clas- 
sical cyclic code [30, 8, 8] with the check polynomial 
h{x) = (1 + + x^), c = 10 and X = 3 gives a code 
[[180, 16, 8]], while X = 1 gives [[180, 16, 6]] with a smaller 
distance. 

Example 12. Same construction starting with the clas- 
sical cyclic code [[30,8,8]] corresponding to the check 



polynomial h{x) = (1 + + x^) with c = 15 and X = 2 
gives a [[120,32,4]] CSS hyperbicycle code; x = 1 gives a 
code [[120,32,2]]. 

Note that in many cases the code rate goes up com- 
pared to the hypergraph-product code constructed from 
the same cyclic codes while the construction from the 
"small" cyclic codes is not possible (cf. Example [5|. 



H. Non-CSS versions of hyperbicycle codes 



We observe that when Hi = Hi and H2 = H2, the con- 



struction in Eqs. (19) can be mapped to non-CSS codes 



in Eq. (16) that in many cases have the same distance 



but half the number of encoded and physical qubits. In 
particular, this happens when x = 1 and matrices Hi and 
H2 are symmetric. By non-CSS hyperbicycle codes we 
then mean a result of the mapping in Theorem [l] of the 
code in Eq. (19). The dimensions of such codes can be 
readily found by applying Theorem |3] where Si = 82 = 0. 

Theorem 10. A quantum non-CSS code constructed 
from matrices (24) such that Hi = Hi and H2 = H2 



and the stabilizer generator matrix 

G = {Ei,0Hi\H2^Ea), 



(40) 



encodes K = Y^ikT'^f'^ /k'i 



€^^^^4^^V4^'^ logical qubits into N = 
cnin2 physical qubits, where pi{x) are all binary factors 

of x^ — 1 such that k^^^ 7^ 0^ including x^ — 1 itself. The 
distance of such a code is bounded by D >\d/c\, d 
mm{di^d2) (same notations as in Theorem 



Proof. The distance bound follows from the proof of The- 
orem [5] given the fact that any code word of the original 
quantum code has to have support on at least one of the 
sublattices with weight exceeding [d/cj. □ 



Theorem 11. Suppose c is even, ai and hi in Eq. I^21y 
are such that /c-"^^^^ = ki, ri = ni and binary codes 
with generator matrices and hi have distance at 

least 2. Then quantum non-CSS code with generators 
in Eq. ( |7^ that have been reduced by construction in 
Eq. IWfhas parameters [[nin2C, A:iA:2, > {2/c)d]] where 
d — min((ii, (i2). 

Proof. This distance bound follows from the proof of The- 
orem [S] given the fact that any code word of the original 
quantum code has to have support on at least one of the 
sublattices with weight exceeding {2/c)d. □ 

Finally, we would like to mention that the upper dis- 
tance bound in Theorem [6] also applies to non-CSS hyper- 
bicycle codes since by construction this bound involves 
only one sublattice. 

For X = 1 we can use palindromic check polynomials 
h{x), i.e. x^®s^'^^^/i(l/x) = h{x), such that cn — d.egh{x) 
is even, in order to construct symmetric circulant matri- 
ces Hi from the polynomial x^'^'^~'^^^^^^^^^'^h{x). 
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Example 13. A [[289,81,5]] non-CSS hyperbicycle code 
(see Fig. ^ is obtained from Eqs. (16) and JT^ using 
circulant matrices Hi = I-L2 corresponding to coefficients 
of a palindromic polynomial h{x) = l-\-x- 
where c = 1 and X = 1 • 



V. CONCLUSIONS 

We described a large family of hyperbicycle codes that 
includes as subclasses the best of the known LDPC codes. 
The construction allows for explicit upper and lower 
bounds on the code distance. We also described new 
LDPC code families with finite rates and distances scal- 
ing as a square root of block length. Our discussion is 
accompanied with geometrical interpretations of the hy- 
perbicycle codes which can facilitate design and applica- 
tions of such codes. The construction is particularly im- 
portant for designing LDPC codes with relatively small 
block lengths which is important since the original hy- 
pergraph product codes have relatively poor parameters 
at small block lengths. 

Another advantage of hyperbicycle construction is that 
it can be based on a pair of very well studied classical 
cyclic codes. This leads to codes with good parameters 
up to limited but relatively large block lengths (in gen- 
eral, cyclic codes with asymptotic rates below one have 
poor asymptotic parameters). The planar layout of thus 
constructed quantum codes possess translational invari- 
ance of stabilizer generators which may simplify the im- 
plementation (see e.g. Ref. [37]). 

Although the quantum LDPC codes discussed in this 
work have been shown to possess a finite noise threshold 



[24| , it is yet to be seen whether there are good decoders 
for such codes. It may well happen that the relation 
between hyperbicycle and bicycle codes can lead to better 
decoding for the former as the latter are known for their 
good decoding properties. 

Even though the lower distance bounds presented in 
this paper are in some cases inferior compared to the 
hypergraph-product codes, we do not expect that this 
will have a significant effect on the value of the noise 
threshold as the distance still scales as a square root of 
the block length while the LDPC structure of the stabi- 
lizer generators is preserved [24 . Given that, we expect 
that one can encode more qubits into hyperbicycle codes 
compared to hypergraph product codes without affecting 
the threshold. 

Our results notwithstanding, there are several open 
questions in regard to the hyperbicycle codes. In partic- 
ular, it would be interesting to establish conditions under 
which the hyperbicycle codes reach the upper distance 
bound. Furthermore, the case when the block shift x 
and the number of blocks c are commensurate has not 
been analyzed. It would also be interesting to explore 
the exact relation between the hyperbicycle codes and 
the CSS codes constructed over higher alphabets [3^ 139]. 
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